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Abstract 



X 

We study the next-to-next-to-leading order (NNLO) evolution of flavour non-singlet quark 
densities and structure functions in massless perturbative QCD. Present information on 
the corresponding three-loop splitting functions is used to derive parametrizations of 
these quantities, including Bjorken-x dependent estimates of their residual uncertainties. 
Compact expressions are also provided for the exactly known, but rather involved two- 
loop coefficient functions. The size of the NNLO corrections and their effect on the 
stability under variations of the renormalization scale are investigated. The residual 
uncertainty of the three-loop splitting functions does not lead to appreciable effects for 
x > 10~ 2 . Inclusion of the NNLO contributions reduces the main theoretical uncertainty 
of a s determinations from non-singlet scaling violations by more than a factor of two. 



1 Introduction 



More than thirty years after the pioneering experiments at SLAC JO, structure functions 
in deep-inelastic lepton-hadron scattering (DIS) remain among the most important probes 
of perturbative QCD and of the partonic structure of hadrons. Indeed, experiments have 
proceeded towards very high accuracy and a greatly extended kinematic coverage during 
the past two decades 0. Moreover, the forthcoming luminosity upgrade of the electron- 
proton collider HERA at DESY will allow for accurate measurements up to very high 
resolution scales Q 2 ~ 10 4 GeV 2 , thus considerably increasing the lever arm for precise 
determination of the scaling violations, i.e., the Q 2 -dependence, of the structure functions. 
An accurate knowledge of the parton densities will also be indispensable for interpreting 
many results at the future Large Hadron Collider at CERN. 

Given the non-perturbative Bjorken-x dependence of the structure functions at one 
scale, the scaling violations can be calculated in the QCD-improved parton model in terms 
of a power expansion in the strong coupling constant a s . The next-to-leading order (NLO) 
ingredients for such analyses are available since 1980 for unpolarized structure functions 
in massless perturbative QCD [0[. Yet the corresponding results for the next-to-next-to- 
leading order (NNLO) are not complete at present, due to the enormous complexity of 
the required loop calculations. Of the components entering the NNLO description, the 
three-loop /^-function (governing the scale dependence of the strong coupling constant) 
H and the two-loop contributions to the coefficient functions (connecting the structure 
functions to the parton densities) || |6], [7|] have been derived. However, only partial results 
have been obtained so far for the three-loop terms of the splitting functions (governing 
the scale-dependence of the quark and gluon densities), most notably the lowest even- 
integer Mellin moments of those combinations relevant to unpolarized electromagnetic 
deep- inelastic scattering ||. 

Standard global analyses of deep-inelastic scattering and related processes, like the 
Drell-Yan process for which two-loop coefficient functions have also been calculated Q, 



have thus been restricted to NLO up to now fllQl , [lT| [Cq| . This level of accuracy is 
however not sufficient to make full use of present and forthcoming data, as the the- 
oretical uncertainties of the NLO results, for instance on the strong coupling constant, 
already now tend to exceed the corresponding experimental errors. Therefore first approx- 



imate NNLO analyses have been performed recently of data on neutrino-nucleon [13] and 



electron (muon)-proton ]TJ[] DIS structure functions, directly using the results of refs. || 
via integer Mellin- N techniques. However, these techniques lack some flexibility, e.g., they 
cannot incorporate additional information on the x-dependence of the two-loop coefficient 



functions [|, 0] and of the three- loop splitting functions [|15|, [16], [17|, |l^, [19[] . Hence we 
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pursue an alternative approach which allows for incorporating the NNLO corrections into 
programs using standard x-space [O, |TTJ or equivalent complex- N techniques . Its 



most important ingredients are compact approximate x-space expressions for the three- 
loop splitting functions including quantitative estimates of their present uncertainty. In 
the present article, we deal with the important flavour non-singlet case. The flavour- 
singlet quantities will be discussed in a subsequent publication. 

This paper is organized as follows: In Sect. 2 we recall the general formalism for the 
scale dependence ('evolution') of non-singlet quark densities and structure functions in 
massless perturbative QCD. The a s -expansions are explicitly given up to NNLO for arbi- 
trary choices of the renormalization and mass-factorization scales. In Sect. 3 we present 
accurate, compact parametrizations of the exactly known || || [7|, but rather involved 
x-dependence of the two- loop coefficient functions. In Sect. 4 we employ the present con- 
straints [||, [T5|, [18] on the three-loop non-singlet splitting functions for deriving approx- 



imate expressions for their x-dependence. The remaining uncertainties are quantified. 
All these results are put together in Sect. 5 to study the impact of the NNLO contribu- 
tions on the evolution of the various non-singlet parton densities and structure functions. 
Here we also discuss the implications on determinations of a s from DIS structure func- 
tions. Finally our findings are summarized in Sect. 6. Mellin-iV space expressions for 
our parametrizations of the two- loop coefficient functions of Sect. 3 can be found in the 
appendix. 
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2 The general formalism 



We set up our notations by recalling the NNLO evolution equations for non-singlet parton 
densities and structure functions. The number distributions of quarks and antiquarks in 
a hadron are denoted by qi(x, /i 2 , /i 2 .) and q~i(x, fij, /i 2 .), respectively, where x represents 
the fraction of the hadron's momentum carried by the parton. yU r and [if stand for the 
renormalization and mass-factorization scales, and the subscript i indicates the flavour of 
the (anti-)quark, with i — 1, . . . , Nf for Nf flavours of effectively massless quarks. 

The scale dependence of non-singlet combinations of these quark densities is governed 
by the (anti-)quark (anti-)quark splitting functions. Suppressing the dependence on x, \i r 
and \if for the moment, the general structure of these functions, constrained by charge 
conjugation invariance and flavour symmetry, is given by 

^yiyfc = ^Mfc = fiikPqq + 'Pqq 

'PqiQk = ^^ilk = ^ikPqq^^Pqq ■ (2-1) 

In an expansion in powers of the strong coupling constant a s the flavour-diagonal ('va- 
lence') quantity V^ q starts at first order, while and the flavour-independent ('sea') 
contributions Vf u , and are of order a 2 . A non- vanishing difference V?„ — V^ ul occurs for 

yy yy ° ^ yy yy 

the first time at third order. This general structure leads to three independently evolving 
types of non-singlet distributions: The evolution of the flavour asymmetries 

QNS,ik = ft ± ft - (ftt ± q~k) (2.2) 
and of linear combinations thereof, hereafter generically denoted by g^s> ls governed by 

^ = K ± K' ( 2 - 3 ) 

The sum of the valence distributions of all flavours, 

N f 

4s = S(ft-ft) , (2.4) 
r=l 

evolves with 

^NS = K - K + N AK - ^S) • (2-5) 

The first moments of V NS and V^ s vanish, 

( 1 dxx N - 1 V^ s = f 1 dx x N ~ 1 V^ s = for N — 1 , (2.6) 
Jo Jo 

since the first moments of g^s an d <?ns reflect conserved additive quantum numbers. 
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The difference V qq — is unknown except for the first moment, which vanishes by 
virtue of Eqs. (|2.3|), ( ^-5|) and (|2.6|). However, the size of the 2-loop contributions to 
and Vg q relative to the corresponding term of V^ q suggests that this difference is negligibly 
small at moderate and large x. Hence we shall use the approximation 

7>ns = ^ns (2-7) 
for the rest of this article, i.e., we henceforth treat q^ s as a '-'-quantity. 

Restoring the dependence on the fractional momentum x and the renormalization and 
mass-factorization scales fi r and /i/, our evolution equations thus read 

r/ , / " 2 

— q^(x, n 2 f , ill) = V^[a s (i4), , 



din/' 2 



2 



Here Cg> stands for the Mellin convolution in the momentum variable, 



[a® 6] (a) = [ ld ^a{y)b(-) . (2.9) 
A y \yj 

The expansion of Pns u P to the third order (NNLO) in a s = otg/ (47r) takes the form 



u 2 

V^ s (x,a s (n 2 r ), -0 = a s (// 2 ) Pis^x) 



-aM (P^W -^(^ln^) (2.10) 



Mr 



+ /5 2 4°s ) (-)ln 2 4)+--- • 

The one- and two- loop functions P^(x) and pjfs( x ) are known for a long time ||; the 
3-loop quantities P^ ± (x) are the subject of Sect. 4. The relevant coefficients of the QCD 



P-f unction, 



are given by 



P(a s ) = -£aj+ a fl , (2.H) 



dln/i 2 



r Z=0 



A = 102 (2.12) 
2857 5033 , r 325 , r2 



& = iVf + Ni . 

H 2 18 } 54 J 



The first two coefficients /?o and (3\ are scheme independent in massless QCD; the result 
given for /?a refers to the MS renormalization scheme employed throughout this paper. 
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The non-singlet structure functions F^ m , a = 1, 2, 3, are in Bjorken-x space obtained 
by convoluting the solution of Eq. ( J2.8D with the corresponding coefficient functions: 



Q 2 m 2 

\lf Mr 



(x) 



(2.13) 



with rji = 2, T] 2 = 1/x, r] 3 = 1, and 

Q 2 i* r 



.->2 ,,2 
A*/ Mr 



5(1 - x) + a s (//^) fc^ s (z) + P§l{x) In 

\ /if 

+« tew + (p^^x) + 



+ 



P (0) (5?) P 
NS ^ ^NS 



(0) 



l''f 

NS ^ c a,NS 
(0) 



(2.14) 



■(°) ^ c (!) 



x) \ In ^- - /?oc£i s (x) In ^ 
J M/ Mr 

2 q 2 ^ 2 " 2 



Q 2 



/iy m 2 



Here an overall electroweak charge factor has been absorbed into g^s- The first-order 
coefficients c^s( x ) can found in ref. ||; the 2-loop quantities c^ 2 J^(x) computed in 
refs. |5], U are discussed in Sect. 3. 

It is often convenient, especially in the non-singlet sector considered here, to express 
the scaling violations of the structure functions in terms of these structure functions 
themselves. The expansion coefficients of the corresponding kernels /C^ NS in 



d 



dlnQ 2 a < NS 



>,Q 2 



/v -a,NS 



Mr 



(2.15) 



are built up of factorization-scheme invariant combinations of the splitting functions 
i"NS ( x ) an< i the coefficient functions c^^si 00 )- Up to third order this expansion reads 



\ Mr / 



(2.16) 



(i) 



C a,NS ^ C a,NS 



-AcSsW+^^ln 2 ^ 

Mr 

- {2/3 {P« ± (x) - ftc^s)} + PiPi 



(0), 



Mr 



This approach removes the dependence of the finite-order predictions on the factorization 
scheme and the scale (if, thus allowing for an easier control of the theoretical uncertainties. 
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3 The 2-loop non-singlet coefficient functions 



The 0(a 2 s ) contributions to the coefficient functions for the structure functions F 2 , 
Fl = F 2 — 1xF\ and xF 3 were calculated some time ago in refs. 0, 0, 0- The resulting 
expressions are rather lengthy and involve higher transcendental functions. Hence it is 
convenient to employ more compact, if approximate, parametrizations of these quantities. 
This holds in particular if one uses the moment-space technique p0"| , which requires the 
analytic continuation of all ingredients to complex Mellin-iV. The reader is referred to 
refs. jn]| for a more rigorous approach to the moment-space expressions for Cf\ Those 
parts of the coefficient functions arising from \x r ^ Q and \ij ^ \i r in Eq. ( |2.14| ) are simple 
convolutions of the well-known lower-order anomalous dimensions and Wilson coefficients. 
The same applies to the terms induced by usual scheme transformations, e.g., that from 
the MS to the DIS factorization scheme. For explicit expressions see refs. 0. Hence the 
parametrizations can be restricted to the MS scheme, and to ln(/i^/Q 2 ) = ln^//^) = 0. 



Our procedure for deriving compact approximate expressions for c a ^ s \ 



We keep the -(—distribution parts, defined by 



dx a(x) + b(x) 



dxa(x) {b(x) — b(l)} 



x) is as follows: 



(3.1) 



exactly (up to a truncation of the numerical coefficients). The integrable x< 1 terms are 
fitted to the exact results for 10~ 6 < x < 1 — 10~ 6 . Finally the coefficients of 5(1 — x) are 
slightly adjusted from their exact values using the lowest integer moments. The resulting 
parametrizations deviate from the exact results by no more than a few permille. This 
holds for the c^ s (x< 1) themselves as well as for the convolutions with typical hadronic 
x-shapes. The adjustment of the 5(1 — x) pieces is important for the latter agreement. 

The non-singlet coefficient function entering the electromagnetic F 2 can be written as 



J2) 



,Ns( X ) 



+ N 



14.2222 L\ - 61.3333 Li - 31.105 L x + 188.64 

1 — x L 

17.19 L\ + 71.08 L\ -660. 7 L x + LiL (-174.8 L x + 95.09 L ) 
2.835 Ll - 17.08 L\ + 5.986 L - 1008 x - 69.59 - 338.046 5(1 - x) 
{ ( 1 



f <; ( : 1.77778 Ll - s, ->926 Ij 4 G.3 189 

- 1.707 L\ + 22.95 L 1 + L!L (3.036 L + 17.97) (3.2) 
+ 2.244 Ll + 5.770 L - 37.91 x - 5.691 + 46.8405 5(1 - x)\ 



with 



L x = ln(l 



x) 



lnx 
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For C2N S (a;), relevant for the charged-current case, the second and third line of this ex- 
pression have to be replaced by 



- 17.19 L\ + 71.08 L\ - 663.0 L x + LiL (-192.4 L x + 80.41 L ) 

— 3.748 Lq - 19.56 L\ - 1.235 L - 1010 x - 84.18 - 337.994 5(1 - x) . 

The corresponding parametrizations for Fl read 



(3.3) 



c L,Ns( x ) 



: 13.62 L\ — 55.79 L\ — 150.5 LiLq + (26.56 x — 0.031)Lq — 14.85 L (3.4) 

+ 97.48 x - 40.41 - 0.164 5(1 - x) H JV> {6 ac£i - 12 xL - 25 x + 6} 

27 



and 



c?,ns(» = 13 - 30 L{- 59. 12 Li- 141. 7 LiL + (23.29 x - 0.043)Ljj - 22.21 L (3.5) 

16 

+ 100.8 x - 52.27 - 0.150 5(1 - x) + — N f {q xL x - 12 xL - 25 x + 6} . 

For the Afy parts, which are identical in Eqs. (3.4) and (3.5), we have taken the exact 
expression from ref. ||. 

Also for the charged-current non-singlet structure function F% there are two combina- 
tions which differ at 0(a 2 s ). The first one, entering F£ N + Fg N , can be written as 



J2)- 
-3,NS 



X 



14.2222 L\ - 61.3333 L\ - 31.105 L x + 188.64 

- 15.20 L\ + 94.61 L\ - 409.6 L x - 147.9 L\L Q 

- 3.922 Lg - 33.31 L 2 - 67.60 L - 576.8 x - 206.1 - 338.635 5(1 - x) 
1 



+ N f l [j—^ [1-77778 L\ - 8.5926 L x + 6.3489 
+ 0.042 L\ - 0.808 L\ + 25.00 Li + 9.684 LiL 



(3.6) 



+ 2.207 + 8.683 L - U.97 x - 6.337 + 46.8575(1 - x) 



For the other combination cj^jt^x), corresponding to Fg N — F£ N , the second and third 
line of this result have to be replaced by 



15.20 L\ + 94.61 L\ - 396.1 L x - 92.43 L\Lq 

3.049 Ll - 30.14 Lg - 79.14 L - 467.2 x - 242.9 - 338.683 5(1 - x) . 



(3.7) 



The complex Mellin moments of these results, c a ^ S (N), can be readily obtained. They 
do not involve special functions beyond the logarithmic derivatives of the T-function. The 
explicit expressions can be found in the appendix. 
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4 The 3-loop non-singlet splitting functions 



Only partial results are presently available for the O(a^) terms P^(x) of the splitting 
functions. In the non-singlet sector, the current information comprises 



(2)- 



the lowest five even-integer moments of P^ s calculated in refs. ||, while for P l 
only the first moment (= in MS) is known; 

the complete 0(NJ) piece (identical for the '+' and '— ' combinations) determined 



via an all-order leading- Nf approach in ref. |T5 



the most singular small- x terms (~ ln 4 x) of P^ + and P^g inferred in ref. 



from the leading small- x resummation of the non-singlet evolution kernels |22j 



The 2-loop results P^s ± ( x ) [0 anc ^ c i 2 Ns( x ) HI furthermore indicate that the difference 
^ns ( x ) — ^ns ( x ) * s negligibly small at large x. Finally present knowledge complies 
with the conjecture that the splitting functions do not receive contributions of the 



form [In (1 — x)/(l — x)]+ with / > in the MS factorization scheme, unlike the coefficient 
functions discussed above. 

In what follows we employ this information for approximate reconstructions of 

Pgt(x) = P o 2)± (x) + N f P? ]± (x) + N 2 f P?\x) . (4.1) 

Our approach is to fix the coefficients of suitably chosen basis functions by the above 
constraints. The spread of the result due to 'reasonable' variations in the choice of those 
functions then provides a measure of the residual uncertainty. Specifically we employ the 
ansatz 

P?W = TT^V + A ^ W ~ x ) + A t* f ^ + A tf^ + A tM x ) + ft (*) (4-2) 
I 1 ~ x )+ 

for the ^-independent (i = 0) and Np (i = 1) terms in Eq. ((4.1[). Here fi and f represent 
contributions which, while being integrable, peak at x —>■ 1 and x —>■ 0, respectively. f m 
stands for a part with a rather flat x-dependence. As for the illustrations in ref. ||, 
these contributions are build up of powers of ln(l — x), x, and lnx. Finally / as allows to 
account for known leading small- X terms. Equating the second to tenth even moments 
of Eq. (|4.2|) to the results of ref. M yields five linear equations which can be solved for 



the coefficients Afj. The case of Pj^ is treated afterwards by taking over A^% and 
from the '-l-'-combinations, as already indicated in Eq. ( [4 . 2| ) , and adjusting the remaining 
coefficients as discussed below. 
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X 

Approximations for the iVj-independent part of P^g , derived from the lowest 
$er moments by means of Eqs. ( }4.2| ) and ( [4.3| ), compared to the exact result. 



Before addressing P^ s (x) we demonstrate our procedure by applying it to a known 
result, the Nf = part P of the NLO splitting function Pns + ( x ) PB"] . In this case the 



leading small- x contributions are lnx and ln 2 x, while the integrable terms most peaked 
at large- x read x 2 and ln(l— x). Disregarding small- x constraints in this example, we thus 
choose 



(4.3) 



The resulting eight approximations are compared to the exact result in Fig. 1 for x < 1. 
The latter curve runs inside the uncertainty band over the full x-range. The moments 
tightly constrain Pq^ + (x) for x > 0.15, the total spread in our approach being about 5% 
at x ~ 0.3. The coefficients of the common leading large-x terms are 



/l( 


x) 


= x 2 


or 


ln(l — x 


fm 


(x) 


= 1 


or 


X 


/o( 


» 


= In x 


or 


In 2 x 


fas 


(x) 


= 







(i> 
o.i 



^0,2 

and the first moments read 



65.13... 68.74 
61. 85... 79.64 



P (1)+ (iV = l) = -2.404. . .0.400 



(exact 66.47) 
(exact 69.00) , 

(exact - 1.127) 



(4.4) 



(4.5) 



'Unreasonable' combinations in the sense of Eq. 



, like ln(l — x), x 2 , and 1 (i.e. 



no 



fo) or 1, lnx, and ln 2 x (i.e., f\ missing), can lead to considerably worse approximations. 
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Figure 2: Approximations for the iVy-independent part of P^g , denoted by P m 
Eq. ( |4.1| ), as derived from the five lowest even- integer moments by means of Eqs. ( [4.2[ ), 
( |4.6| ) and (|4.7p . The full lines represent those functions selected for further consideration. 



(2) + 

Now we turn to P ( x )- The additional loop or emission may, besides adding two 
powers of lnx, lead to two additional large-x logarithms with respect to Pq 1 ^ + (x) (the 
transition from 1-loop to 2-loop yields however only a term ln 1 (l — x) ). Hence we put 

fi(x) = x 2 or ln(l — x) or ln 2 (l — x) or ln 3 (l — x) 

f m (x) = 1 or x 

fo(x) — In x or In 2 a; (4-6) 

/ as (x) = ^C|(ln 4 x + Aln 3 x) . 

Besides A = we also include A = —4 and A = 8 for / = hi a;. Subleading small- x 
terms of this order of magnitude are suggested by the expansion of Pns + an d Pns + i n 
moment-space around = [25]. Thus we consider 32 combinations, 8 of which are 



rejected as they fail to fulfill the further ad hoc, but mild constraint 

100P (1)+ (iV = l) < P (2)+ (iV = l) < -40P (1)+ (iV=l) (4.7) 

on the perturbative expansion of the first moment. The x> 1 behaviour of the remaining 
24 function is displayed in Fig. 2; their 1/(1 — x) + coefficients span the range 

1138 < A { q} < 1625 (1347) . (4.8) 

The bracketed number applies if combinations with fi(x) = ln 3 (l — x) are disregarded. 
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Due to the larger function pool of Eq. ( |4.6| ), the large-x uncertainty band of Fig. 2 
is some factor of three wider than that for P (a;) in Fig. 1, reaching a total spread 
of about 15% at x ~ 0.3. Moreover Pn 2 (x) is rather unconstrained at small x ^ 10~ 2 



by present information, as the leading small- x term [|T8[ does not dominate over less 
singular contributions at practically relevant values of x. However, physical quantities are 
only affected by the splitting functions via convolutions with smooth non-perturbative 
initial distributions which 'wash out' the oscillating large- x differences of Fig. 2 to a large 
extent. Furthermore the convolutions receive important contributions from the (well- 
constrained) large- x region of Pns(^) even at very small x. The above 'bare' uncertainty 
is thus considerably reduced over the full x-range. This effect is illustrated in Fig. 3, 

(2)+ 

where four representative approximate results for P are convoluted with a simple, 
but typical input shape. The total spread after this convolution is as small as 0.3% for 
0.2 < x < 0.9, and becomes large only at x < 0.02. 

The uncertainty band of Fig. 3 is rather completely covered by the results 'A' and 'B'. 

(2)4- 

Hence our final estimates for P ( x ) an d its remaining uncertainty are given by 
P^(x) = 1137.897 ^ + 1099.754 5(1 - x) - 2975.371 x 1 



x 



125.243 - 64.105 hT x + 1.580 In 4 x (4.9) 
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P (2 l + (x) = 1347.207- — + 2283.0115(1 - x) - 722.137 ln 2 (l - x) 

(l-x)+ 

- 1236.264- 332.254 In x + 1.580 (In 4 x- 4 In 3 x) . (4.10) 
The average \ [ Pq 2 ^~ (x) + Pq 2 ^ (x) ] represents our central result. 

The iVj-term P^ + is the leading radiative correction to P 1 ( ' 1 ' )+ (x), which is in turn 
only slightly more complicated than the 1-loop non-singlet splitting function p3fl . Hence 



it is natural to adopt here the ansatz ( |4.3|) employed for the 2-loop ./V°-piece in our above 

(2) + 

illustration. The resulting eight approximations for P{ (x) are displayed for x < 1 in the 
left part of Fig. 4 (dashed curves). Their spread at large x is similar to that obtained for 
Pq^ + in Fig. 1. The leading large- x coefficients fall into the range 

- 190 < A {2 \ < -180 . (4.11) 

(2)+ 

The uncertainty of the complete result for P^s ( x ) * s dominated by the spread of the 
above iV/-independent contribution, as estimated by the difference between Eqs. (fO|) and 
( |4.10| ). This is also true at small x, despite the fact that the band in Fig. 4 is presumably 
an underestimate in this region, as a possible term ~ In 3 x has been disregarded. Hence 
it is sufficient, at the present stage, to keep only the iV°-contribution to the error band 

(2}+ (2)+ 

of Pjsj S and to employ just one representative for P x . Our choice, an average of two 
typical results with and without a ln(l — x) term, reads 

p} 2)+ (x) = - 184.4098 - — 180.6971 6(1 - x) - 98.5885 ln(l - x) 

(1 - x) + 

+ 205.7690 x 2 + 6.1618 + 5.0439 In 2 x (4.12) 
and is also shown in the left part of Fig. 4 (solid curve). 



As mentioned before the iV 2 -piece in Eq. (f4.1|) is exactly known from ref. After 
transformation to x-space, this contribution reads 

P?\x) = 1 ( - 64 - [204 + 192 C(3) - 320 C(2)] 6(1 - x) + 64 

+ (96 In x + 320) + (1 - x) (48 In 2 x + 352 In x + 384)) , (4. 13) 

1 — x J 

where ((I) denotes Riemann's (- function. 

Finally we consider P^~(x). Here our treatment is inevitably more approximate. 
According to the expectations given at the beginning of this section, we take over the 
1/(1 — x) + and 6(1 — x) terms of the '+ '-combinations in Eqs. (|4.9|), ( |4.10| ), and ( [4.12| ). 
The remaining coefficients are (after inserting the appropriate leading small-x piece 
||18||) determined by the first, eighth and tenth moments of ref. [EL assuming that the 
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Figure 4: Left: approximations to the Nj part of P$s + i obtained from the five lowest 
even- integer moments using Eqs. (|4.2| ) and Q4.3] ). Right: approximate results for the N® 

and Nj terms of the 3-loop splitting function P r 



(2)- 
NS • 



difference to P^ s (x) is negligible for the latter two, entirely large-x dominated quantities. 
The results are shown in the right half of Fig. 4. The uncertainty band for P (x) is 

(2) + 

about 50% wider than that for P ( x ) around x = 0.3, reflecting the lack of precise 
knowledge of the intermediate- N moments, but smaller at small-x, as this region plays a 
much greater role for the first moment known here from Eq. (2.6), than for the second 
moment in the '+'-case. Our parametrizations spanning the present uncertainty are given 
by 

1 



1137.897- 



+ 1099.754 5(1 -x)- 2954.684a;' 

2 



1 — X t 

- 143.709 - 2.761 In 2 x + 1.432 In 4 x 
1347.207- — - 



x 



+ 2283.011 5(1 -x)- 722.238 ln 2 (l - x) 



(4.14) 



(4.15) 



- 1234.756 - 327.479 In x + 1.432 (In 4 x - 4 In 3 x) 
(also here the average represents the central result), supplemented by 

Pi (2) "(x) = - 184.4098 - — l — 180.6971 5(1 - x) - 98.5722 ln(l - x) 

(1 - x)+ 

+ 205.3670 x 2 + 6.5740 + 3.5474 In 2 x . (4.16) 

(2)4- 

For the latter expression an average has been calculated in the same manner as for Pf J . 
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5 Numerical results 



We are now ready to consider the numerical impact of the NNLO terms on the evolution 
of the non-singlet parton densities and structure functions. Before doing so, however, it 
is worthwhile to look at the perturbative running of a s underlying these considerations. 
In the left part of Fig. 5 the a s -expansion (2.11) of the /3-function is shown for Nf = 4 
flavours. Besides the contributions of Eq. (2.12) relevant for NNLO calculations, also 
the contribution — @zol b s of ref. p6| has been included. If one uses the effect of this four- 



loop (N 3 LO) term as an estimate of the residual error of the expansion, the resulting 
uncertainty amounts to 0.08%, 0.35%, 1.1% and 2.5% for a s = 0.12, 0.20, 0.30 and 
0.40, respectively. The effects are somewhat larger (smaller) for Nf = 3 (Nf = 5). The 
consequences of this expansion on the scale dependence of a s are illustrated in the right 
part of Fig. 5. For this illustration we have used Eq. (2.11) with Nf = 3 at fi r < m c = 1.5 
GeV, Nf — 4 between m c and = 4.5 GeV, and Nf = 5 for fi r > nib, assuming that 
a s {f^r) i s continuous at these thresholds. If a s is fixed to 0.115 at fi r = Mz, then the four- 
loop effect reaches 0.1% (1%) only at /i 2 = 20 GeV 2 (1.5 GeV 2 ), respectively. Clearly the 
truncation of the series (2.11) after three terms does not introduce a significant theoretical 
uncertainty in the kinematic regime of deep-inelastic scattering. 
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Figure 5: Left: The perturbative expansion of the QCD /5-function up to order a^, for 
four flavours in the MS renormalization scheme. Right: Illustration of the resulting scale 
dependence of a s , using a variable Nf as detailed in the text, /z 2 is given in GeV 2 . 
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For illustrations of the scale dependence of the parton densities and structure functions, 
initial distributions have to be chosen at some reference scales, in the following denoted 
by and Qq, respectively, in Eqs. (2.8) and (2.15). We will employ the function 

/ = x°- 5 (l-x) 3 (5.1) 

for all six quantities 

/ = xq£ s (x,fj, 2 ffi ), F± NS (x,Q 2 Q ) and xF± m (x,Q 2 ) . (5.2) 

Eq. (5.1) represents a simple model shape which incorporates the most important features 
of non-singlet ^-distributions of nucleons. The same input is used in all cases, as this allows 
for a direct comparison of the effects of the various kernels in Eqs. (2.10) and (2.16). The 
overall normalization of / is irrelevant for the logarithmic derivatives considered below. 
Our initial scales are specified via 

a«(^r=^/,o) = a s(v 2 = Ql) = 0.2 , (5.3) 

irrespective of the order of the expansion. For a s (M|) = 0.114. . .0.120 this choice cor- 
responds to — Qq ~ 25 ... 50 GeV 2 , a Q 2 -region typical for fixed-target DIS. If not 
explicitly indicated otherwise, the results will be given for Nf — 4 massless flavours. 




Figure 6: The perturbative expansion of the scale derivative, q£ s = d\nq^ s /d\nii 2 , for a 
non-singlet '-(-'-combination of quark densities at fi r = fif. The initial conditions are as 
specified in Eqs. (5.1)-(5.3). Here and in what follows the subscripts A and B indicate 
the approximations for the 3-loop splitting functions derived in the previous section. 
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The evolution of q^ s (x,fij) is illustrated in Fig. 6 for the standard choice /i r = [if 
of the renormalization scale. In this case the perturbative expansion appears to be very 
well convergent: Except for the region around x ~ 0.07 where the scale derivative is 
very small, the NNLO corrections for q^ s = d\nq^ s /d\n fj.'j are as small as about 2%, 
while the NLO contributions typically amount to 10 . . . 20%. The residual uncertainty 
of the 3-loop splitting functions of Sect. 4 leads to a noticeable effect only for x < 0.02, 
and even at x ~ 10" 3 this effect does not exceed ±2% with respect to the central result 
^(NNLOyi + NNLOs) not shown in the figure. Over the full x-range the NNLO corrections 
are comparable to the dependence on the number of active flavours: If Nf is increased 
(decreased) to Nf = 5 (Nf = 3), q^ s is decreased (increased) by about 2%, respectively. 
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Figure 7: The dependence of the NLO and NNLO predictions for dlnq^s/dlrL/if on the 
renormalization scale ji r for six typical values of x. 
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Another way to assess the reliability of perturbative calculations is to investigate the 
stability of the results under variations of the renormalization scale fj, r . In Fig. 7 the 
consequences of varying ji r over the rather wide range < fif. < 8/ij are displayed 
for six representative values of x. The relative scale uncertainties of the average results, 
estimated by 

= max [g£ s (x, $ = \n) ■ ■ ■ - min [q£ s (x, = \n) . . . 
qm ~ ' ~ 2 | average [g+ s (x, ^ = \fi . . . 4^)] | [bA) 

are shown in the left part of Fig. 8. Also this estimate leads to about 2% for the NNLO 
uncertainty, an improvement by more than a factor of three with respect to the cor- 
responding NLO result. Even as low as x ~ 1CT 3 the NNLO calculation, despite its 
approximation uncertainty increasing towards small x, is superior to the NLO. 

Finally the evolution of ' — '-combinations g^s * s illustrated in the right part of Fig. 8. 
For x > 0.1 the difference to the '+'-case discussed so far is negligible at NLO as well 
as at NNLO. At small x the NLO predictions differ by up to 2%. As expected from the 
discussion in Sect. 4, the residual uncertainty of the NNLO result is considerably less 
pronounced at small x in the '—'-case, but somewhat larger for 0.01 <x < 0.1. 
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Figure 8: Left: The renormalization scale uncertainty of the NLO and NNLO predictions 
for the scale derivative of q^ s , as obtained from the quantity Ag^s defined in Eq. (5.4). 
Right: The NNLO effects on the evolution of q-^ s for the standard scale choice \i r = /if, 
together with a comparison of the NLO partonic NS + and NS~ evolutions. 
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We now turn to the evolution of the non-singlet structure functions. The physical scale 
derivative F 2 ~ns = dlnF^s/dhiQ 2 is shown in the left part of Fig. 9 for /i r = Q. Besides 
the splitting functions Pj$s ( x ) the effect of which has been illustrated in Fig. 6, here also 
the coefficient functions c^lss{ x ) an d c^^(x) enter the NLO and NNLO evolution kernels 
as detailed in Eq. (2.16). These additional terms considerably increase the Q 2 -dependence 
at large be seen by comparing Fig. 6 and Fig. 9. E.g., the NNLO corrections rise 

from 4% at x — 0.5 to about 7, 11 and 21% at x — 0.65, 0.8 and x = 0.95, respectively. 
The corresponding NLO contributions amount to 24, 30, 37 and 51% of the LO results. 
Unlike for the parton densities, the NNLO corrections to the structure functions are larger 
than the ^/-dependence at large x: If Nf is increased (decreased) to Nf = 5 (Nf = 3), 
P 2 n S is decreased (increased) between 3.5% and 7% for 0.5 < x < 0.95, respectively. 

The worse convergence of the expansion at large x is due to the large soft-gluon 



contributions [ln fc (l — x)/(l— x)] + , k 



21-1, to c 2 



(0± 

NS 



x: 



which are conjectured to 



be absent [24] in the MS splitting functions P } 



an 

NS 



(X 



Consequently, as shown in the right 
part of Fig. 9, keeping only the coefficient-function contributions in the O(a^) term of 

/•o\_i_ 

Eq. (2.16) yields a very good approximation at large x. In fact P^ s contributes less than 
2% to the total NNLO derivative P 2 "ns at x > 0.2. The residual uncertainty of Pj^c, + \ 
given by the difference NNLO a - NNLO#, is thus completely negligible in this region. 
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Figure 9: The perturbative expansion of the scale derivative, F 2 ns = dlnF^s/dlnQ 2 , 
for a non-singlet structure function at // r = /if. The initial conditions are as specified in 
Eqs. (5.1)-(5.3), the leading-order curve is identical to that of Fig. 6. Also shown (right 

( 2) + 

part) is the effect of omitting the contribution from the 3-loop splitting function P^g . 
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The dependence of -FVns on the renormalization scale /i r is presented in Fig. 10 and 
Fig. 11 (left part), analogously to the partonic case (see Eq. (5.4)) in Figs. 7 and 8 using 



AR 



2,NS 



max [F+ NS (x, = \Q 2 . . . 4Q 2 )] - min [F 2 + Ng (x, = \Q 2 . . . 4Q 2 
2 | average [F+ NS (x, = \Q 2 . . . 4Q 2 )] | 



. (5.5) 



The slower large-x convergence of the a s series for F 2 "ns * s obvious from these results as 
well, e.g., no extremum close to /j, r — Q is obtained for x = 0.8. The NNLO uncertainties 
as estimated using Eq. (5.5) read 3%, 4.5% and 7% for x = 0.5, 0.65 and 0.8. The 
corresponding NLO results are 8.5%, 10.5% and 12%, respectively. The accuracy of the 
Q 2 -slope predictions is thus improved by a factor 2 ... 3 except for very large x. 
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Figure 10: The dependence of the NLO and NNLO predictions for dlnR^g/dlnQ 2 on 
the renormalization scale ji r for six typical values of x. 
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Figure 11: Left: The /^-uncertainty of the scale derivative of -F 2 + ns; as estimated by 
Ai^fs defined in Eq. (5.5). Note that the absolute values of F^ s are ver y sma U f° r 
0.05 < x < 0.15. Right: The NNLO effects on the evolution of -F 2 7ns f° r Mr = A*/, together 
with a comparison of the NS + and NS~ evolutions for -F 2)NS at NLO. 




Figure 12: The scale derivative -F 3 7ns = d\nF 3 ~ NS /d\nQ 2 at /i r = /i/. The initial condi- 
tions are as given in Eqs. (5.1)-(5.3), the leading-order curve is the same as in Figs. 6 
and 9. Also shown (right part) is the ratio of F 3 ~ m to the corresponding result for F^ NS . 
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As for the parton densities shown in Fig. 8, the evolution of -F 2 ~ns illustrated in the 
right part of Fig. 11 is indistinguishable from that of -F 2 "ns a ^ x > 0.1, while being better 
constrained at NNLO at very small x. For x ~ 10~ 3 , the 2.5% . . . 3.5% positive effect of 

(2) — 

Pjsjg in Fig. 8 is overcompensated by the coefficient-function contributions. This effect 
also occurs for F^ s n °t displayed at small x. In both cases the NLO corrections are 
smaller than for q^ s , resulting in a better small-a: NLO renormalization-scale stability of 
F^ NS as can be seen by comparing the left parts of Fig. 11 and Fig. 8. 

The scaling violations of F£~ NS are presented in Fig. 12 for the medium- to large- a; 
region. Since the soft-gluon terms [ln fc (l — x)/(l — x)] + are identical in c^jJs and cjjjJs, 
the results for F 3; ns and F 2 ^s agree (for identical initial distributions as assumed here) 
as x — > 1. However, the different regular terms lead to noticeable differences already 
at medium x, reaching 5% and 10% at x ~ 0.4 and 0.3, respectively. At small x the 
corrections are considerably larger for F^s than for F2,ns> resulting in scale uncertainties 
of about 10% at NLO and 4% at NNLO for 10~ 3 < x < 10~ 2 for the former quantity. 

Finally we turn to the determination of a s from scaling violations of non-singlet struc- 
ture functions. Here we address the uncertainties Aa s which arise from the truncation of 
the perturbation series, confining ourselves to the region x > 0.25 of considerable negative 
scale derivatives i^.ivs- I n this region the results for i^Ns are rather similar to those for 
^2,ns and need not to be considered separately. We also disregard the negligible large-x 
differences between the scaling violations of -F 2 "ns and F 2 ~ m and between the NNLOyi 
and NNLO# calculations. Our procedure for estimating Aa s is as follows: For each x 
we determine those scales [i r ,m\n and ^ r>max which led to the minimal and maximal NLO 
and NNLO results for |^2,ns| use d in Eq. (5.5). The value of a s (Ql) is then adjusted to 
obtain, at these values of x and /i r , the same results for -F 2jNS as found for \i r = Q and 
tts(Qo) = 0-2 (Fig- 9, left part). The latter standard-scale results thus play the role of the 
experimental results for F 2i ns m determinations of Aa s in data fits. 

The resulting upper and lower limits for a s {Q^) are shown in the left part in Fig. 13. 
Due to the increase of the higher-order corrections towards large x discussed above, the un- 
certainty Aa s rises with increasing x. As available experimental DIS results are restricted 
to x < 0.85 pi, we choose a value x ~ 0.55 for estimating the x-averaged uncertainties 
given by the differences to the reference result a s (Ql) = 0.2. This procedure yields 

Aa s (Q2~ 25...50 GeV 2 ) NLO = > (5-6) 

A« s (g 2 ~ 25...50 GeV 2 ) NNLO = ^OtM ■ ( 5 - 7 ) 

Often results and uncertainties for a s from different processes and observables are com- 
pared after evolution to a common reference scale, conventionally chosen as the Z-boson 
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mass Mz- Adopting Ql = 30 GeV 2 (and Nf = 5 for Q < \i r < Mz) for definiteness, one 
obtains the error bands displayed in the right part of Fig. 13 and 



0.006 



0.0025 



Aa s (M|) NLO - Iq"qq4 > Aa s (M|) NNLO - I £0015 ■ ( 5 - 



As expected from our previous discussions below Eq. (5.5), the NNLO calculation reduces 
the theoretical uncertainty under consideration by a factor of about 2.5. 

In a data analysis, also the NLO and NNLO central values for a s (Ql) will be different, 
since the NNLO scaling violations are stronger over most of the large-x region as shown 
in Fig. 9. A simple estimate analogous to that for Aa s yields 



a s (M|) NNLO -a s {M 2 z \ 



NLO 



0.002 



(5.9) 



Due to the strong x-dependence of the NNLO/NLO ratio, this estimate is less reliable 
than Eq. (5.8), its uncertainty amounts to about ±0.001. Nevertheless it is interesting 



to note that Eq. (5.9) agrees with the findings of ref. [I3| from analyses of data on F 3 . 
The 3-loop splitting function P^s contribute only about —0.0007 to the shift (5.9) of the 
NNLO result. 
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Figure 13: Right: The x-dependent theoretical uncertainty of the determination of a s 
from the scale derivative of -Fj.ns at Ql = 25 ... 50 GeV 2 , estimated by the fi r - variation 
\Ql < A*r ' — 4Qo- The scales leading to a maximal (minimal) |F2Ns| are denoted by /i rimax 
(// r ,min)- Left: The resulting error band for a s (M§) using Ql = 30 GeV 2 and Nf = 5. 
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6 Summary 



We have investigated the effect of the NNLO perturbative QCD corrections on the scale 
dependence of flavour non-singlet quark densities and structure functions. For this pur- 
pose, and for application in further analyses, we have derived compact parametrizations 
of the corresponding three-loop splitting functions P^g and the two-loop coefficient func- 
tions Ca^s, a = 2,3, L. The latter quantities are exactly known [g, ||, [JJ; their analytic 
x-dependent expressions are however rather cumbersome and not readily transformed to 
moment space [21\. Our parametrizations of cjj^ s (x) and their Mellin transforms thus 
provide a convenient technical tool. They agree to the exact results up to a few permille 
or less over the full x-range, thus introducing a negligible error of well below 0.1% after 
insertion into the perturbative expansions. 

As only partial results are presently available for the three-loop splitting functions 
H |15|, our parametrizations of P^g (x) serve the additional purpose of providing 
quantitative estimates of their x-dependent residual uncertainties. The function P^ s (x), 
relevant to the evolution of flavour asymmetries like u + u — [d + d] , is well constrained 
at large x by the lowest even-integer moments of refs. H, the spread reaching about 
±7% at x ~ 0.3. On the other hand PjSjg (x) is very weakly constrained for x < 10~ 2 so 
far, despite the fact that the leading small-x term is known [|T^]. The quantity P^g (x), 
entering the evolution of the quark-antiquark differences, is somewhat better (worse) 

(2)4- 

constrained at small x (medium x), respectively, than P^ s (x). As the splitting functions 
enter parton densities and structure functions only via convolutions with smooth non- 
perturbative initial distributions, these 'bare' uncertainties are very much reduced for 

(2)4. 

physical quantities over the whole x-range. E.g., the spread of P^g leads to effects of 
less than ±0.2% at x > 0.2 after convolution with typical nucleonic input shapes. In this 

(2)4- 

region the present uncertainties of P^g are thus rendered absolutely negligible, leading 
to effects even below 0.01% after insertion into the perturbation series. Their impact 
becomes significant only for x < 10 -2 , without seriously impairing the NNLO calculations 
even down to x ~ 10 -3 . 

The perturbative expansion for the scale dependence dlng^x, /i 2 )/<i In /z 2 of the non- 
singlet combinations of quark densities appears to be very well convergent. For a s = 0.2, 
corresponding to scales of about 25-50 GeV 2 , the NNLO effects of Pjj^ are on the level 
of 2% rather uniformly in x. This result is to be compared to the NLO corrections which 
amount to 10-20%. Also the variation of the renormalization scale leads to effects of 
about ±2% at NNLO. Corrections of this size are comparable to the dependence of the 
predictions on the number of quark flavours, rendering a proper treatment of charm effects 
rather important even for large-x non-singlet quantities. 
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Especially at x > 0.5, the higher-order corrections are much larger for the scale deriva- 
tive dlnF aNS (x,Q 2 )/d\riQ 2 , a = 2,3, of the non-singlet structure functions. This en- 
hancement is an effect of the coefficient functions containing large [hi fc (l — x)/(l — x)] + 



soft-gluon terms, which are conjectured to be absent in the MS splitting functions |21 . 
E.g., the NLO and NNLO effects reach about 37% and 11% of the respective lower-order 
results at x = 0.8 for a s = 0.2 and four flavours. The NNLO calculations thus represent a 
distinct improvement, reducing also the renormalization-scale dependence of the predic- 
tions by a factor of two to three, e.g., to about ±7% at x = 0.8. Accordingly the inclusion 
of the NNLO corrections into fits of data on non-singlet scaling violations is expected 
to yield, besides a slight lowering of the central values for a s (M|) by roughly 0.002, a 
considerable reduction of the (so far dominant) theoretical error due to the truncation of 
the perturbation series, 

A +0.0025 
Aa s (M z J N NLO - _ 0.0015 

These estimates are compatible with the results of the fits of F 3i Ns-data performed in 



ref. 13 , where an alternative, integer-moment based approach to the calculation of the 



scaling violations has been pursued. 

Fortran subroutines of our parametrizations of c a ^ s (x), a = 2,3, L, and P^ s (x) 
can by obtained via email to neerven@lorentz.leidenuniv.nl or avogt@lorentz.leidenuniv.nl. 
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Appendix: Third-order quantities in Mellin- AT space 



The Mellin transforms of the approximate NNLO expressions of Sect. 3 and Sect. 4 are 
given in terms of the integer- N sums Si(N) and their analytic continuations 

N 1 f-lV 
Si=S l (N) = J2- l =ai)-±-^^\N+l) . (A.l) 

Here £(1) stands for the Euler-Mascheroni constant, and £(Z > 1) for Riemann's £-func- 
tion. The I th logarithmic derivative i/>( i-1 ' of the T-function can be readily evaluated 
using the asymptotic expansion for Re N > 10 together with the functional equation. 

Due to the simplicity of our parametrizations for the three-loop splitting functions, 
only the most simple Mellin transforms occur for these quantities. Therefore we are able 
to dispense with details here. The Mellin- N dependence of the exactly known iV^-piece 



can be found in ref. 



The moments the non-singlet '+ '-coefficient function (3.2) entering F 2 e ' m ' are given by 

c?ns ( N ) = + 3-55555 (6 S 4 + 8 S 3 Sx + 3 Sf + 6 S 2 Sl + Sf 

+ 20.4444 (2 S 3 + 3 S 2 S l + Sf) - 15.5525 (S 2 + Sf ) - 188.64 S 1 
165.356 /15.38 9.7467\ 358.503 2.9678 3 

+ "iv - 53 " {-jf —) S2 + ^^ S2Sl + ^r Si 

/ 174.8 9.7467N 2 /190.18 116.734N 

+ [ipr + ^r) s i - {—z —) s ^ 

17.01 34.16 306.849 72.5824 1008 nnnnii 

H — H = 338.044 (A. 2 

N 4 N 3 N 2 N N+l y ' 

+ N f {- 0.59259 (2 S3 + 3 S 2 Si + Sf) - 4.2963 (S 2 + S 2 ) - 6.3489 Si 

6.072 /6.072 18.0408\ /6.072 17.97 14.3574\ 

~ nr Ss " vw ~n ~ vw + "^v - ) Sl 

0.07078 o2 4.488 4.21808 21.6028 37.91 inninn 

+ ^ s l + ^r + —! jv FTT + 468406 

For the charged current '—'-combination the third to fifth line of this result are, according 
to Eq. (3.3), replaced by 

229.916^ /31.58 9.7467\ „ 393.703 „ „ 2.9678 

/192.4 9.7467X 2 /160.82 61.1321N 

+ 1^5-+ IC ' 



22.488 39.12 265.774 164.777 1010 
+ ^ iV^ + ^ iV iV+T - 337.992 . 

The first two lines and the sixth line of Eq. (A. 2) stem from the universal -(—distribution 
parts of Eq. (3.2). They are exact up to a truncation of the numerical factors M. 
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The corresponding A^-space results for the coefficient functions (3.4) and (3.5) for Fl read 
(2)+ 136.88 13.62 c2 ^55.79 150.5N 0.062 

(A ° " " "Iv -52 + Hv -51 + w) Sl "W 

14.85 207.153 53.12 97.48 A /A . 

+ IF- + — + (FTTP + FTT " 0464 (A - 4) 

16 f 6 _6_ __6 25_ 

+ f 27 I ivTT + jv + (at + 1) 2 ~ ]v+T 



(2) _ 128.4 13.30 2 / 59.12 141.7\ 0.086 
cim{N) = -^-S2 + ^-S 1 + [— ^ Inl- 



and 



N N \ N N 2 J N 3 

22.21 180.818 46.58 100.8 „„„„ 

H =-H hi r^H 0.150 (A.5) 

A 2 (N + l) 3 N+l V ; 

16 f 6 6 6 25 

+ '27 1 ?+i 1 + n + (ivTTp ~ ivTT 

Here the A 7 / parts represent an exact result [0. 

The ' — '-coefficient function (3.6) for F3, occurring in the u+v sum, leads to 

4 2 ns(^) = + 3.55555 (6 S 4 + 8 S 3 Si + 3 Sf + 6 S 2 S 2 + Sf) 

+ 20.4444 (2 S 3 + 3 S 2 £i + 5?) - 15.5525 (S 2 + S?) - 188.64 S 1 
297.756 „ /147.9 33.2767\ „ 298.733 „ „ 0.9778 ^ 

( 147.9 33.2767\ 2 45.8683 

23.532 66.62 67.6 373.029 576.8 /A , 

+ — — + — — -FTT" 338 - 625 < A - 6 > 

+ iV> I — 0.59259 (2/33 + 3 S 2 Si + Sf) - 4.2963 (s 2 + Sf) - 6.3489 ft 
°-° 42 (2 5s + 3 ft ft + Sf) + °^S! + f ^ - ft 



— ^ ^ + 6 + + ^r-^r + ^ — ; 

10.6538^ 4.414 8.683 15.9177 14.97 

+ N~ 2 + ~Jp W N FTT + 4 

For the '+ '-combination of Eq. (3.7) entering F% N — F 3 uN one has to replace the third to 
fifth line of the above result by 

186.816 „ Z92.43 33.2767\ „ 187.793 „ „ 0.9778 ^ 
Z92.43 33.2767\ o2 123.121 



+ — j^ft (A.7) 



+ { N 2 + N 

18.294 60.28 79.14 276.473 467.2 

+ — — + -Tp ~ — ~ NTT ~ 338 ' 681 
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